Geodesic stability for KS Black hole in Horava-Lifshitz gravity 

via Lyapunov exponents 
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By computing the Lyapunov exponent, which is the inverse of the instabihty time 
scale associated with this geodesic motion we show that there is two region of space 
which in both of them the equatorial timelike geodesies are stable via Lyapunov 
measure of stability. 
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MATHEMATICAL PRELIMINARIES:LYAPUNOV EXPONENTS 

Let us consider a trajectory (in phase space) described by a certain evolution. The 
Lyapunov exponents (also known as characteristic exponents) associated with a trajectory 
are essentially a measure of the average rates of expansion and contraction of trajectories 
surrounding it. They are asymptotic quantities, defined locally in state space, and describe 
the exponential rate at which a perturbation to a trajectory of a system grows or decays 
with time at a certain location in the state space [1]. Analysis conducted with Lyapunov 
exponents are called Lyapunov stability analysis. They are useful in characterizing the 
asymptotic state of an evolution (attractors in dissipative systems) [2]. Using Lyapunov 
exponents, we can distinguish among fixed points, periodic motions, quasiperiodic motions, 
and chaotic motions [3]. 

Concept of Lyapunov Exponents 

We begin by defining Lyapunov exponents for a given system of equations. Let X{t) 
such that X(t = 0) = Xq represent a trajectory of the system governed by the following 
n-dimensional autonomous system: 

x = F{x-M) (1) 

where the vector x is made up of n state variables, the vector function F describes the 
nonlinear evolution of the system, and M represents a vector of control parameters. Denoting 
the perturbation provided to X{t) by y{t) and assuming it to be small, we obtain an equation 
after linearization in the disturbance terms. The perturbation is governed by 



y = Ay 



(2) 
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where, in general, A = D^F[x{t); M] is an n x n matrix with time dependent coefficients. If 
we consider an initial deviation y{0), its evolution is described by 

y{t) = 0(i)y(O) (3) 

where is the fundamental (transition) matrix solution of (3) associated with the tra- 
jectory X{t).The eigenvalues of A provide information about the stability of the associated 
fixed point. 

The procedure used to determine Lyapunov exponents can be considered to be a generaliza- 
tion of linear stability analyses. An interesting and detailed discussion on the relationship 
between linear stability analysis and Lyapunov stability analysis can be found in the paper 
of Goldhirsch, Sulem, and Orszag [4]. They argue that the Lyapunov exponents are global 
quantities associated For an appropriately chosen y{0) in (3), the rate of exponential expan- 
sion or contraction in the direction of y{0) on the trajectory passing through Xq is given 

by 

A,- = limt^^ J log (4) 

where the symbol 1 1 denotes a vector norm . The asymptotic quantity Xj is called the 
Lyapunov exponent. We have n Lyapunov exponents associated with an n-dimensional 
autonomous system. It can be shown that if the trajectory X{t) corresponds to a motion 
other than a fixed point, then one of the Xj is always zero . Following Lyapunov [5], the 
fundamental matrix (f){t) is called regular if 

limt^^j log det{(f){t)) (5) 

exists and is finite. If 0(t) is regular, then, according to a theorem , the asymptotic quantity 
defined in (4) exists and is finite for any initial deviation y(0) belonging to the n-dimensional 
space. 



STABILITY OF CIRCULAR ORBITS VIA LYAPUNOV EXPONENTS 

As was shown by Cardoso et al[6 ] for all spherically symmetric spacetimes, in a geometri- 
cal optics approximation, QNMs can be interpreted as particles trapped at unstable circular 
null geodesies and slowly leaking out. The leaking time scale is given by the principal Lya- 
punov exponent, which one can obtain a fairly simple expression, in terms of the second 
derivative of the effective radial potential for geodesic motion. As was stated by Cardoso et 
al, this deep, intuitive approach to the QNMs and it's relation to the circular null geodesies 
is valid for all asymptotically fiat, spherically symmetric black-hole spacetimes. The most 
important note is that the powerfuU formalism which is presented in [5] is completely inde- 
pendent from the form of the action, i.e there is no difference between treating a stationary, 
asymptotically flat, static spacetime from a higher order gravity as f(R) or HL and the usual 
higher dimensional (or 4-dim) BHs in GR.Thus with have no worrying about the validity of 
their method we can use it and following all the steps of the [6]. Without no loss of generality 
we can restrict ourselves to a simple kind of problems includes circular orbits in stationary 
spherically symmetric spacetimes and equatorial circular orbits in stationary spacetim .This 
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case could be described by a two-dimensional phase space which for our geodesies analysis 
may be described by 

X{t)^{pr,r) (6) 

Linearizing the equations of motion about orbits of constant r , the principal Lyapunov 
exponents can be expressed as 

Where Vr — P and a dot denotes the derivative with respect to the proper time r. Remember 
to mind that For circular geodesies 



KS BLACK HOLE SOLUTION IN HL THEORY 

There are different versions of the HL theory. As a geometrical point of view all these 
versions Follow from the ADM decomposition of the metric [7], the fundamental objects of 
interest are the fields N{t,x), Ni{t,x), gij{t,x) corresponding to the lapse , shift and spatial 
metric of the ADM decomposition. In the (3 + 1) -dimensional ADM formalism, where the 
metric can be written as 

ds"^ = -N^dt"^ + Qijidx' + N'dt){dx^ + N^dt) 

and for a spacelike hypersurface with a fixed time, its extrinsic curvature K^j is 

Kij = ^{g\j-V^N,-VjN,) 

where a dot denotes a derivative with respect to t and covariant derivatives defined with 
respect to the spatial metric gij, the action of Hofava-Lifshitz theory for 2; = 3 is 

[ dtd^x^N{CK- Cv) 

we define the space-covariant derivative on a covector t>j as ViVj = diVj — F'^v; where F'^- 
is the spatial Christoffel symbol, g is the determinant of the 3- metric and N — N{t) is a 
dimensionless homogeneous gauge field. The kinetic term is 

Here Ni is a gauge field with scaling dimension [iVj] = 2. 

The potential term Cy of the (3 -|- l)-dimensional theory is determined by the principle of 
detailed balance [8] , requiring jCy to follow, in a precise way, from the gradient flow generated 
by a 3-dimensional action Wg. This principle was applied to gravity with the result that 
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the number of possible terms in Ly are drastically reduced with respect to the broad choice 
available in an 'potential is 



4A — 1 



Where in it Cij is the Cotton tensor [9 ] which is defined as, 

The kinetic term could be rewritten in terms of the de Witt metric as: 

2 



K 



hi 



Where we have introduced the de Witt metric 

Gijkl^]^i^gikgjl^gilgjk^_y^gijgkl 

Inspired by methods used in quantum critical systems and non equilibrium critical phe- 
nomena, Hof ava restricts the large class of possible potentials using the principle of detailed 
balance outlined above. This requires that the potential term takes the form 



^2 

— -^E'^^GijkiE' 



Note that by constructing E^^ as a functional derivative it automatically transverse within 
the foliation slice, ViE^^ = 0. The equations of motion were obtained in [10]. KS BH is a 
static spherically symmetric solution for HL theory which contains 2 parameter , one mass 
like parameter m and a parameter which controls the escape from a naked singularity ou and 
satisfies [11] 

2 1 
cum > - 

- 2 

In the usual spherical coordinates {t, r, 9, (p) and in the Schwarzschild's gauge the metric 
reads: 

ds^^diag{+f,-j,-r^j:2) (8) 
where in it the metric gauge function is 



/ = 1 + our'^ - \JuJ^r'^ + Amur (9) 

and is S2 the surface element on a unit 2- sphere. As motivated by Sekiwa "it is obvious 
that l/2ui is equivalent to and this means that we could view l/2uj as a charge in some 
degree" [12]. Thus The outer and inner event horizon can be compared with the outer and 
inner event horizon of Reissner-Nordstrom black hole [13]. Essentially as claimed by the 
founders of the KS, this solution "represents the analog of the Schwarzschild solution of 
GR\ 
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CIRCULAR ORBITS 

The best treatment of the geodesies equations is due to Chandrasekhar [14]. For a simple 
typical form of a spherically symmetric metric in Schwarzschild gauge 

ds^ = diag(+f(r), -r^Es) (10) 

one can obtain the following expression for the potential function: 

yr = 9ir){^^-^-6,) (11) 

Where in it the E, L respectably can be interpreted as the Energy, angular momentum of 
the test particle within the circular orbits in similar to the classical mechanics and 6i = 1,0 
for timelike and null geodesies, respectively. For KS BH (8) the relation (11) convert to: 

Vr^E'- - SJir) (12) 

TIMELIKE GEODESICS 

For circular orbits we know that both the potential term and the first derivates of it must 
be vanish which leads us to the next expression for the second derivative of the potential 

thus using (7) the Lyapunov exponent at the circular timelike geodesies is 

A = ^V(2/-r/0K" (14) 



Since the energy must be real, we require 

aiog(/(r)) 



< 2 (15) 



d\og{r) 

this inequality leads to the next bound for the radius of the timelike circular orbit ^ 

oo^r^ - 9m^uj^r + Amuj > (16) 

This is a cubic algebraic equation . We review the general solution for a cubic eq.Let the 
general cubic eq be 

x^+px + q^O (17) 
To have three distinct real root (the case which is happen in KS BH) we define 

A = V + 27q^ (18) 



^ Remember that in KS BH we must have m^w > 5. 
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If A < then these three roots can be obtained by solving a simple trigeometric 

cos(3^) = cos(q;) (19) 



where in it 



for radial eq (16) we have 



where in it 



cos(a) = (20) 
X — A cos{9) 



A = -108{-yd < -297(-)2 (21) 

OJ (jO 



5 = 27m^uj^ - 4m^c^^ > 2.75 (22) 
finally we can write the next expressions for all roots of the (16) 

r„ = 2V3mcos(— -±-),n-0,l,2 (23) 

Also we know that the KS solution has a physically accessible horizon located at 

r = /i = m + y/m^ - (24) 

The KS solution is vahd only for such region of the space where r > h. We must determine 
that which of these three distinct real values r^, i = 0, 1, 2 belongs to this interval .Indeed we 
must have 

h<rn< 2VSm (25) 

(Which is satisfied automatically) 

Also the negativity of the energy impose that for timelike geodesies we must restricted 
ourselves to 

(ro < r < ri) Y (r > r2) (26) 

This bound for the radial coordinates which arisen from the reality of the Energy has an 
essential role for investigating the stability via the relation (14) for the Lyapunov exponent. If 
we want to have a stable motion the exponent A must be imaginary and have a negative 
real part and unstablee unless. For instability we must have 

V;' > (27) 

The dominator of it must be checked.That is we must check that when r belongs to the (26) 
the next inequality satisfied or not? 

-3ff/r + 2f"-ff>0 (28) 

the opposite sign indicates on astable one. A very simple and straightforward calculation 
show that for all values of the r belong to the (26) we have 

K" < (29) 
Thus the time like geodesies in KS BH are stable. 
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CONCLUSION 

HL theory brings some important new features from the GR to the higher dimensional 
lagrangian and it's role in construction a non relativistic candidate for quantum gravity.The 
Lyapunov exponents is a very important and powerful method to stability any dynamical 
system both in classical mechanics and also in quantum treatment of complex systems.! 
n this short letter after a brief review of the Lyapunov method in mathematical physics, 
following the method which have been described by Cardoso et al[5] for 4 and 5 dimensional 
BHs in GR,we applied a similar method to the a new type of the BH which indeed are 
a formal IR limit of the HL theory and show that the timelike geodesies are stable under 
small perturbations. As Cardoso et al show that for a equatorial circular timelike geodesies 
in a Myers-Perry black-hole background are unstable .But we show that in the context of 
the HL theory these geodesies are stable. Thus we can argued that the nonlinear terms are 
successful to pass from the instabilities. 
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